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Abstract 

We consider methods of interpolating between the crystalline and dynamically triangulated random sur- 
face models. We argue that actions based on the deviation from six of the coordination number at a site are 
inadequate and propose an alternative based on Alexander moves. Two simplified models, one of which has 
a phase transition and the other of which does not, are discussed. 
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Many apparently disparate physical systems, such as non-critical string theory, 2D quantum gravity, 
and polymer membranes, are described by the statistical mechanics of random surfaces (for recent reviews 
see [1]). This letter is concerned with understanding the differences, if any, in the critical behaviour of 
crystalline surfaces and dynamically triangulated surfaces. 

Crystalline surfaces are described by a fixed triangulation T (in continuum terms a fixed intrinsic 
metric), with the partition function formed by summing over all the possible embeddings of this surface in 
a D-dimcnsional euclidean space [2]; usually the fixed triangulation is a regular one but this is not essential. 
If the lattice sites i of T have coordinates Xj in the Z?-dimensional embedding space the action is given by 

St{k,D) = ± Yl (X l -X,) 2 + K ]T 1-Aa.Aa' (1) 

<ij>GT AA'eT 

where < ij > denotes the link fron i to j and a' are the un it normal vectors of the triangles on cither 
side of the link. The canonical partition function is 



Z T ( K , D) = TT / d D X, S D VxJ e-^(^) (2) 




where the delta function suppresses the translational zero mode. The first term in the action is the elastic 
contribution while the second term which is proportional to the extrinsic curvature of the surface inhibits 
bending. The model exhibits the crumpling transition, characterised by the Hausdorff dimension of the 
random surface changing from dn = 2 (the smooth phase) to dn — > oo (the crumpled phase) at some critical 
value of the extrinsic curvature coupling. The existence of this transition, which appears to be of second 
order, is not proven analytically but is well established by the numerical studies [2,3,4,5]. 

In the case of a dynamically triangulated surface (DTRS), the canonical partition function is formed 
by summing over all possible embeddings and all possible distinct triangulations T with a fixed number of 
points N keeping one point marked 

Z n (k,D)= ]T P (T)Z t (k,D) (3) 

T:\T\=N 

where p(T) is a weighting factor depending only upon the local properties of T [6,7,8]. The grand canonical 
partition function is then 

oo 

Z( f i,K,D)=J2 e ~' jN M^D) (4) 

JV=1 

The sum over triangulations is equivalent to the integral over worldsheet metrics in Polyakov's ([9]) formu- 
lation of string theory; hence these models are used in the numerical study of non-critical strings and 2-D 
quantum gravity. It is known that in order to obtain an interesting continuum limit in which both the string 
tension and the mass gap scale to zero it is necessary (but not necessarily sufficient) to have the extrinsic 
curvature term in the action [10,11]. In the last few years much computational effort has been expended 
on these models but the results of these simulations are much less clear-cut than for the crystalline sur- 
face. The Monte-Carlo simulation of systems with dynamical triangulation is very computationally intensive 
which limits the size of lattice that can be studied and the difficulties are compounded by the fact that the 
crumpling transition, if it exists, seems to be rather weak with only a cusp singularity in the specific heat 
[12,13]. 

As an alternative to studying the DTRS directly we consider a generalized model which interpolates 
between it and the better understood crystalline surface. Introduce the notion of a distance, n(T ,T), 
between a general triangulation T and reference triangulation T (we shall consider shortly what form this 
distance might take) and consider the canonical partition function for a generalized model 

Z' n (C,k,D)= p(T)exp(-(n(T ,T))Z T (K,D) (5) 

T:\T\=N 
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When C = this coincides with the DTRS while as £ — > oo all triangulations except To are suppressed and 
we obtain the canonical partition function for the crystalline surface on To. It is slightly more complicated to 
define a grand canonical partition function for this new model because not all forms of reference triangulation 
T) can be written naturally for all N. However, in the cases where a natural definition is possible, such as 
the example we will discuss later, the grand canonical partition function is given by 

oo 

Z'( f j l ,C,K,D)=Y,e-'* r Z' N (<;,K,D) 

N=l 

At large ./V the asymptotic behaviour of Zn is expected to be 

Z N = AH" + O (J-j^J (7) 

Accordingly the grand canonical partition function develops non-analytical behaviour as \i J, [i c and the 
critical exponent 7 controls the derivative of Z in this limit through 

dZ 1 

The phase diagram of the generalized model in the (k, () plane certainly depends on D; in practice, 
most work on the crystalline surface and the DTRS has been done in D = 3. It is known that the crystalline 
surface has a crumpling transition at k — k c ; for n < k c a surface of N sites has mean square radius 
R 2 ~ log A whereas for k > k c , R 2 ~ N. Since there is no smooth interpolating function between ./V and 
logiV when N — > 00 these two regions must be separated by a phase transition. The simplest possibility 
is that there is a line of transitions (the solid line in fig.l) extending all the way across the diagram and 
cutting the ( = axis at k' c ; in this case the DTRS also has a crumpling transition. However there may 
be a critical point at X so that the order of the transition is different for the DTRS; there is no numerical 
evidence for a first order transition but a higher order one is certainly possible. Alternatively, the line may 
stop at an intermediate point X provided that one or both of the lines of transitions AB, CD is present to 
isolate the two phases of the crystalline surface (the order of A, C in the picture is of no significance and they 
might coincide with X); in this case there is no crumpling transition in the DTRS but there may be some 
non-trivial behaviour arising from the proximity of X to the £ = line (it was pointed out in [13] that the 
observed behaviour is quite similar to that seen in models which are known not to have a phase transition 
but rather display some crossover behaviour). Furthermore, since it is certainly known that R 2 ~ N v with 
v w when (k, Q — (0,0), if the crumpling transition line does stop at X then the line of transitions CD 
must exist. Assuming that actually v = 0, the existence of the line AB cannot be deduced from the mean 
square radius; however there are other quantities, such as 7, which can differentiate between phases. The 
rest of this paper is concerned with the behaviour along the k = line. 

A recent paper [14] proposed using the coordination number of the vertices as a measure of the dis- 
tance n(T, Tq). In a regular triangulation all (or almost all, depending upon the genus) the vertices have 
coordination number <fc = 6; thus we could choose 

n(T,T ) = 5>i-6| a (9) 

for some positive power a to give the distance of T from a regular triangulation. In [14] it was argued that 
a > 2 corresponds to an irrelevant operator in the continuum limit and so a = 1 was chosen. However, it is 
clear that for a > 2, n(T,T ) increases faster as T becomes more different from T than it does for a = 1; 
thus it is hard to see how a = 1 can be more effective at suppressing the irregular triangulations regardless 
of considerations about the continuum limit. Let us assume that the suppression of irregular triangulations 
is strong enough that when £ — > 00 only the most regular possible triangulations contribute. For genus zero 
there are many triangulations with {qi — 6, Vi} differing only by their modular parameter which can vary 
from t = 1 (ie the shortest cycle is O(VN) links long) to r = O(N) (ie the shortest cycle is 3 links long); 
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(6) 



(8) 



they all contribute to Zn and are not suppressed at all by the C,n{T, To) term in the action. Thus even at 
D = the action (9) is not suitable to interpolate between the fluid and the crystalline surfaces. At non-zero 
D the X integral in (2) can be done when n = to get [6,7,8] 

Z N (0,D)= ]T p(T)(Dct7 T )~* (10) 

T:\T\=N 

where Tr is the incidence matrix of T. At large N we know that 

(Det'/^o^)))" 1 ~ e? N (Det'/^D)" 1 (11) 

for finite positive (3 [15]. This means that in the thermodynamic limit the ultra-thin configurations will 
dominate in (10). Hence, even in the limit £ — > 0, the partition function will be dominated not by a nice 
regular two-dimensional surface but rather by quasi one-dimensional spaghetti-like objects. We believe that 
this phenomenon may be responsible, in part, for the observation in [14] that, although the intrinsic geometry 
is smoothed by (9), the extrinsic geometry is not. Clearly the long thin tubes will be suppressed by the 
extrinsic curvature term when n > 0. However, the extrinsic curvature can certainly not be relied upon to 
suppress all the triangulations except one; in addition the k region where the spaghetti freezes out could get 
mixed up with the crumpling transition of the r w 1 triangulations and add to the confusion. 

To make progress a better definition of n(T, T ) is needed. It is known that any two triangulations T , T 
with the same number of vertices and the same topology are related by a sequence of link flip operations 
(known as Alexander moves). Let us define n(T, T ) to be the minimum number of flips necessary to change 
T into T . (It is doubtful that the absolute minimum is really needed; the number of flips taken by any well- 
defined reproducible procedure will do as well.) The new definition of n(T, To) will get rid of the spaghetti. 
It takes O(N) flips to turn a r = 1 triangulation into one with t = O(N) so the suppression factor e~^™ 
can always be made to overcome the effect of (11) by choosing £ large enough. This definition of n(T, To) 
seems much harder to use than (9), especially in a numerical simulation; however there are some special 
cases (unfortunately in D = 0) which can be solved and we now turn to them. 

It is convenient to work in terms of the dual graphs rather than the triangulations themselves. We first 
consider three point functions and take as our reference graph Gq(N) the one shown in fig. 2 (note that the 
actual number of points in the graph is 2A*" + 1). It was demonstrated in [16] that any three point function 
graph can be reduced to this form by flips. Of course such a reference graph is very different from the one on 
which we would like to define the crystalline surface but it serves to illustrate the principle. In the sum over 
the graphs we include all those which can be constructed from tree graphs in the manner shown in fig. 3; Tfc 
denotes the set of rooted trees with k branches and it is known [16] that the number of graphs, T k , in this 
set is given by the fc'th Catalan number 

T (2fc-2)! 

It is easy to see that the rooted tree T^ in a graph of the form shown in fig. 3 can be reduced to the form 
of the reference graph by making k — 1 flips, each time flipping the link that is the trunk of a tree. Now 
introduce the generating function 

"max 

F N {t) = J2 n N ,nt n (13) 
n=0 

where fijv.n is the number of graphs which can be moved to the reference graph Gq(N) by making n flips 
and t = e~^. The graphs obey the Schwinger-Dyson equation shown in fig. 4 and hence Tjv(i) obeys the 
recursion relation 



N 

\ ^ _,fe- i - 
k=l 



F N (t) = yj*- 1 T k F N _ k (t) (14) 
Multiplying by z N (z = e _At ) and summing over N gives, after some manipulation, 

oo oo 

g(z, t) = J2 F *( t ) = 1 + *) 5> f ) fcTfc ( 15 ) 



N=0 fc=l 
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The last sum is just the generating function for the Catalan numbers so we find that the grand canonical 
partition function 9(z,t) is given by 



t < 
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Zc = 


1-t 
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1 




Zc = 
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Zc = 




2' 




4i 



x 2* - 1 - VI - 4zt 
G ™ = 2(z-l + t) (16) 

^(x, i) displays non-analyticity as we take z — > z c (i) to obtain the thermodynamic limit; we find the following 
behaviour 

1 

z c = 1-t 

G~(zc-z)-i (17) 

dg , ,1 

The behaviour for t > \ corresponds to 7 = \ which is the value expected for a system of branched polymer- 
like configurations; given the method of construction of our set of graphs that is not surprising. For i < 5, 
Q itself diverges; at t = it is easy to check that we get exactly the partition function expected when only 
the reference graphs Gq(N) are included demonstrating that varying t does indeed allow us to interpolate 
between an ensemble of graphs and a fixed reference graph. The standard thermodynamic quantities are 
easy to calculate and we find 

l(„( G , Go) ) = {A |*j m 

and 

±MG,G«f)c = {f? «{ (19) 

There is no divergent specific heat or other behaviour to suggest that a correlation length is diverging at the 
critical point t = 5. 

The model can be generalized slightly to include all graphs made of 7^ as shown in fig. 5 (considering 
the four point function in this case simplifies the calculation). By making a succession of flips as indicated 
this can be turned into a graph in which the root of Tk is isolated as before. The number of flips to do this 
is 2min(<7 — l,k — q) so that the factor Tk in (14) is replaced by 

k 

T k X ^^min^-l.fc-g) (20) 
q=l 

(In fact this does not give us the minimum number of flips to convert these graphs to G form but it does 
yield a soluble model.) Similar manipulations as before lead to the result 

g(z t) = t 2 (l-t 2 ) ^^^^^^ 

{,) t 2 (l -t 2 )-t (VI - 4zt 2 - VI - 4it) + |(1 - tf (VI + 4zt 2 - VI - 4zt 2 ) 

Despite its apparent complexity the non-analytic behaviour of Q is the same for all t and takes the form 

Q{z,t) — (22) 

Z c — Z 

There is no phase transition; there is no need for one because 7, the only order parameter that could 
distinguish between t = and t = 1 actually takes the same value in each case. 

Although our definition of n(T,T ) has some appealing properties it is hard to see how it might be 
implemented efficiently in a Monte Carlo simulation to study the evolution of the crystalline surface into the 
dynamically triangulated one. To take a random triangulation and "flatten" it into a regular To by means 
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of flips in the minimum possible number of moves seems to be a problem of the travelling salesman type. It 
is easy enough to count flips going away from To so long as their concentration is small but this is hardly 
likely to be the interesting regime for the transition between the two surface models. 
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Figure Captions 

1) Illustrating possible phase diagrams for the generalized model. 

2) The reference dual graph Gq(N). 

3) Restricted set of three point graphs constructed from tree graphs. 

4) Schwinger-Dyson equation for the restricted set of three point graphs. 

5) A larger set of graphs constructed from tree graphs. 
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